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Qh' Abstract 



a 



Generalizing the linearisation procedure used by Dirac and later by Levy-Leblond, we 
jrt ' derive the first-order non-relativistic wave equations for particles of spin 1 and spin 3/2 starting 

'~{ ' from the Schrodinger equation. 

^^, By the introduction in the momentum of a correction linear in coordinates, we establish 

the wave equation of the radial harmonic oscillator with spin-orbit coupling. 

1 Introduction 

fSJ , The standard procedure to derive the Dirac's free particle equation T] is based on the search of 

^^ ' linear relation between momentum and energy starting from Klein-Gordon equation. Making use 

r — , of the Dirac idea, Levy-Leblond [2] performed such linearization procedure for the Schrodinger 

f^ ' equation and established a linear equation in both momentum and energy which leads to the Pauli 

equation with the correct value of the Lande g factor. A generalisation of such equation was realised 
for higer spin particles using the Bargmann-Wigner method in Galilean Relativity and the Galilei- 
invariant wave equation [3] , [3] and [S] . In [5] Hurley has proposed, in a rigorous and sophisticated 
manner, a Galilean higher spin wave equation. Using the Galilei group and from general invariance 
assumptions, the author has established a first-order wave equation with (6s + 1) components for 
C^ ' spin-s particles which admits a consistent quantum-mechanical interpretation. 

In this work, we are interested by the question of why the Dirac and Levy-Leblond linearization 
algebraic procedure leads only to spin 1/2 wave equation. Is it possible to generalise this method 
to obtain wave equations for particles of spin greater than 1/2. 

On the other hand, we focus on the non-relativistic harmonic oscillator. The Schrodinger's 
equation for this system involves a quadratic potential in coordinates and it's seems natural in the 
context of first-order equations to search a linear potential describing this system. Such problem 
has been studied in relativistic quantum mechanics. The Dirac oscillator system [6]- [9] leads to 
the standard harmonic oscillator in the non-relativistic limit, the Duffin-Kemmer-Petiau oscillator 
[ID]- [13 is a generalisation to spin zero and one particles. This task has been realised in the 
framework of the five dimensional Galilean covariance in [18] and |19j . 



X 



Our purpose in this work is of double interest: the first is the generahsation of the Dirac and 
Levy-Leblond procedure to spin 1 and spin 3/2 particles. The obtained equations are not new 
but coincide with those established by Hurley. In the present paper, the derivation of the wave 
equation is achieved without referring to Galilean invariance and group theoretical techniques and 
it's based on the linearization concept. The second goal is to introduce in these wave equations a 
linear potential in coordinate which leads to the standard non-relativistic harmonic oscillator plus 
the spin-orbit coupling for spin 1/2, 1 and 3/2 particles. 

2 Linearisation 

In the Schrodinger equation, the presence of the energy operator E is linear while the momen- 
tum p appears quadratically in the Hamiltonian. To eliminate this asymmetry between the time 
(d/dt) and spaces variables {d/dx), we require that the wave equation must be of first order in all 
derivatives. This leads to put the wave equation in the form 

{AE + B.p + C)ij = (1) 

here A , B {Bi, B2, B3) et C is a set of five operators with no dependence on E and p and which 
have to be determined. If we act on the left of ([T]) with the operator 

A'E - B.p + C (2) 

where A' and C" are a two new numerical matrices, we find 

{A'AE^ + {A'B, - B,A)p,E + {A'C + C'A) E - B^B.p^, + {B,C - C'B,)p, + C'C) tA = (3) 

the repeated indices are summed over i,j = 1,2,3. If we now impose on ([3]) to be identical to 
Schrodinger's equation 

{2mE - p2) ^ ^ (4) 

the following relations must be satisfied 

A' A = C'C = 0, A'B, - B,A = 0, (5) 

A'C + C'A = 2m, C'Bi - BiC = 0, (6) 

B,B,p,p,=p\ (7) 

Let us note that there is 7 unknown operators to be found A , Bi, C, A' and C unlike ten operators 
in the case of Levy-Leblond procedure. We will see in the next sections that the expressions and 
dimensions of these numerical matrices are related to the spin of the particle. 



3 Particle of spin 1/2 



We start by solving the last set of algebraic equation in the spin 1/2 case and we will reproduce 
the Levy-Leblond equation by a way which could be generalised to spin greater than one-half. The 
ai Pauli matrices obey the identity 

CTiajPiPj = {pi) l2x2, (8) 



where l2x2 is a 2x2 unit matrix. Comparing ([7|) and (|H]), it is straightforward to show that 
Bj ^ aj, because otherwise the second equation in ([5]) would lead to A' = A which contradicts 
the first conditions in ([5]). We look for operators Bi depending on Ci solutions of ([7]). Since in ([5]) 
there is a product of two cr^ matrices with different indices i and j, we choose the form 

which leads to ([7]). The determination of Bi makes the task easier for the rest of the operators. 
Indeed, using ([5]) and dH), one obtains 

here 1 is the unit matrix. Taking into account © , we see that is sufhcient to put A + ^ equal to 
1 . Combined to the first equation in ([S]) , this leads in one hand to A' + A = -^ + ^ = 1 ^^'^ o^ 
the other hand to A = ^ and A' ^ 4^. The following solutions are then obtained 

The wave function ij] — \ I is now a four components object where f and x are two-component 

spinors. Finally, the first order wave equation is written as 

i:y-{ii)''Mi'i)V-'- "" 

which is nothing else but the Levy-Leblond equation reproduced here in a different way. The 
Clifford algebra doesn't appear as in the original paper of this author. However, we will see in the 
next sections devoted to the spin 1 and 3/2 particles that this method can be easily generalized. 

Now, we will introduce in (jl3|) the harmonic oscillator potential. We seek some expression that 
should be linear in r and leading to the three-dimensional harmonic oscillator. We choose the 
substitution 

p ^>- p — imLurjr, (14) 

where w denote the frequency of the oscillator, r the position vector (r = Xi) and r; is a matrix 
equal to 2A^ — 1. The motion equation for this system becomes 

{AE + B.{p-imLO-qr)+C)i:^Q, (15) 

which turns into a coupled system of equations for the components Lp and x 

2mx = —(T.{p~imujr)ip (16) 

Eip = —(T.{p + imujr)x, (17) 



multiplying P^ hy a {p + imuir) and using (fT7)) and the identity 

aiaj ~ Sij + ieijkCFk 
we obtain the wave equation for the Lp component 



Elp 



where 



P 1 2 2 3 ^ 2w , ^ , 

2m 2 2 ?i ^ ^ 



r X p. 



<P, 



(18) 



(19) 



Equation (118^ corresponds to the standard harmonic oscillator with the addition of the spin-orbit 
coupling. This result is in agreement with that obtained by the non relativistic limit of the Dirac 
oscillator [Bj. 



4 Particle of Spin 1 

The problem is to construct the five unknown operators A , Bi and C obeying conditions ([S]) to 
([7]) for a such particle. Following an approach similar to the one used above for spin 1/2, we start 
by solving ([7]). The Bi operators in this case, should contain the 3x3 spin one matrices which are 
given in the standard representation by 



{Si)]k = -«feyfc 



(20) 



£ijk being the totally antisymmetric Levi-Cevita symbol. Inserting the identity —ih^^ i^-P) u=p x 
u where u — {ui,U2, U3) is a vector field, into the well-known vector identity V (V.m) — V x (V x 
u) ~Au leads to a relation between the square of the momentum and the three components of the 
spin operator, i.e. 



h ^ {N-jN^ + SjSi) PjPi = p^l3x3 

where N^ denote the transposed matrix of Ni which is a 3 x 3 matrix defined by 



N,, 



J3xl 



JSxl 



with 



ei = (1 0), 62 = (0 1 0),e = (0 1). 
In the same way as in the spin 1/2 case, one may write for the Bi matrices 



B, 




Thus, from 1^^ one obtains 
BjBiPjp^ 



h-'^ 








OjOi 






PjPi 



(21) 

(22) 
(23) 

(24) 



(25) 



The non diagonal terms are zero 

{Nfp,) {s,pi) = {s,p,) {N,pi) = (26) 

but the diagonal elements are different from zero and consequently equation ^ is not satisfied. 
This is due to the presence in (|2ip of an additional term contrary to equation ([5]) for the spin 1/2 
case. Imposing in the identification procedure that the spinor ip must be included, the condition 
([7]) becomes after acting on ip as 

B,B,p,p,^ = pV- (27) 

In the present case, the wave function i/) is a three components object given by 






-0 = V with M = U2 \ ^v — \ V2 and ti; = W2 \ ■ (28) 
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It is straightforward to show that (|27| is satisfied when 

p.w = (29) 

t;=t> (i'i,W2 = 0,W3 = 0) . (30) 

This equations means that the third component of ip has a zero divergence and that the second must 
be a scalar. We will see that these conditions will be automatically fulfilled when the remaining 
operators A and C will be determined. Using (P^ and in an analogous way as in the spin one half, 
one can easily obtain the numerical expressions of these matrices. Thus, the first-order free wave 
equation for spin one particle takes the form 

srs\ /000\" 

= 0, (31) 

where N— {Ni,N2,N3). Separating the components, one immediately verifies the relations (P5)) 
and ([50]) . Since the matrices N^ have two rows composed entirely of zeros, ([?T|) can be reduced 
to an equation with 6s + 1 = 7 components and then coincides with the wave equation established 
by Hurley 5j using a different approach centered on the Galilean invariance. 

We now examine the harmonic oscillator problem. We introduce this potential via the following 
substitution 

p — )■ p — imuj-qr, (32) 

where the matrix rj is given now by 

/ 1 \ 
77 = 2^2-1 = -1 . (33) 

V -1 ; 

Equation (|3ip can be decomposed into a set of the coupled equations 

Eu — —hr^ [N {p + iniujr) v + s {p + imcur) w] 

2mv = —h^^N^ (p — imur) u (34) 

2mw = —h^^s {p — imojr) u, 



thus the wave equation for the field u is given by 

2mEu = h-^ [(^•P+) (N^-P-) + {s-P+) (s-P-)] u (35) 

where p^ — pzLimcur. From the definitions (pn|) and 1^^, one readily shows that Ni and Si verify 



N,N^ 



S^Sj — tflO'ij^S}^ 



h'5. 



With the help of the above relation, the evaluation of (|55|) yields 



£;« 



P , -^ 2 2 

— TOW r 

2to 2 



-/iw 



CJ 



(L.s) 



(36) 



(37) 



which describes the usual isotropic harmonic oscillator plus the spin-orbit interaction with the 
strength {—uj/H). Let us note that this strength is one half that obtained for spin 1/2 case and 
it's coincides with the non relativistic limit of the DKP oscillator [TU] and that obtained in the 
framework of the five dimensional Galilean covarience ^9] . 

5 Spin 3/2 particle 

For such particle, we will adopt the following matrix representation for the components of spin 
operator s 



si 
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We seek the Bi solutions of ^ containing the hermitian matrices Si. We introduce the 4x2- 
matrices Ki and their 2x4 hermitian conjugates Kf defined by [20] 



Kt 



,K+ 











^3^ = Ho 72 oj 



The Ki and Si matrices and the components of the the momentum p obey the relation 

(4/9) h-^ {KjK+ + Sjs,) p0, = p^Uy^i. 



(39) 



(40) 



Since the previous equation is similar to (I?!]) , one can take the form given by ([M)) where Ni and h 
are replaced by Ki and 3h/2. As in the spin one case, we find that the (6s + 1) x (6s + 1) following 
matrix 

KjK++SyS, 

BjB,p,p,^{A/9)n-^ \ K+K, K+s, \ pjp„ (41) 



3 ' 

„2 



is not equal to p^ and therefore ([7]) is not satisfied. Again, the requirement that the wave function 
must be taken into account leads to ()27p . where now f/; is a ten-component object 

f 
V- = I ^ I , (42) 

X 



here (/? and x are four-component and i7 is a two-component functions. Equation (|27p yields for 
the spin 3/2 particle 

(4/9) h-' [{K,p,) (K+p,) + {s,p,) (s,p,)] ^ = pV, (43) 

(4/9) ft-2 (i^+p^.) [(;^^p^) f7 + (5^p^) ^] ^ p2j^^ (44) 

(4/9) h-^ {s,p,) [{K,p,) n + (s^p,) x] = p'x. (45) 

The first equation is obvious and the last two equations leads to conditions on the components of 
the wave function 

n^~ (2a/3) {K+p,) ^ (46) 

X = - (2a/3) (s,K) ¥>, (47) 

where a is a constant to be determined. The last step is to construct the A and C matrices 
consistent with the two previous equations. We proceed with a method different to the one used in 
the previous cases and without using ([5]) to ([7]) deduced from the linearization procedure. Taking 
into account the expression of Bi and (jH]), the resolution of equation 

B,p,,l,^^{AE + C)^, (48) 

subject to the constraints (|46| . (j47|) and (|4]), we find that A and C are given by expressions 
similar to that of spin one case with the appropriates dimensions and the a constant is given by 
a — (2raK)~^ . Finally, the (6s -I- 1) components first-order wave equation for this particle is given 

by 

100\ I Q K s\ /OOO 

£; + (3n/2)"M K+ p + 2m 1 

000/ \ s {) {) I loot 



i'^Q. (49) 



This result coincides with the Hurley's equation |5] for spin 3/2 particle. 

For the harmonic oscillator problem, again we perform the substitution p —^ p ~ imujrjr with 
rj = 2A^ — 1. The decomposition of the wave aquation into its components and the elimination of 
VL and X i^i favour of Lp gives 

2m^(^ = (4/9) n-2 [(K.p+) {K+.p_) + (s.p^) (s.p_)] ^ (50) 



with the help of the property [20] 

K,K+ + s,s^ = I (3ft/2) e,jkSk + (9/i/4) S,, (51) 

obtained using the definitions (|38l) and (p9]) , the evaluation of the right-hand side of (|50|) yields 



Eifi 



P^ 1 2 2 3^ 2w , ^ , 

2to 2 2 3h^ ' 



V. (52) 



which corresponds to the isotropic harmonic oscillator plus the spin-orbit coupling whose strength 
is one third to the one obtained for spin 1/2 case. 

6 Conclusion 

A new method is presented based on the linearization of the wave equation which allow us to 
deduce the first order equations for non-relativistic particles of spin 1/2, 1 and 3/2. Requiring that 
the wave function ip must be taken into account in the set of relations ([5]) to (O, we were able 
to construct the operators A , Bi, C for particles of spin 1 and spin 3/2. This requirement is no 
longer necessary for a particle of spin 1/2 and the resolution of the previous set of equations leads 
then to the particular solution established by Levy-Leblond. 

In addition, we considered the problem of the harmonic oscillator in the framework of linear 
wave equations and we were able to find the standard wave equations plus the terms giving the 
spin-orbit coupling in each case. 

The method presented in this paper can be generalized to the case of a non-relativistic (and 
relativistic) particle of arbitrary spin as well as that of a spinless particle (works in progress). 
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